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Abstract: We develop a new method to study the thermalization of time depen¬ 
dent retarded Green function in conformal held theories holographically dual to thin 
shell AdS Vaidya space times. The method relies on using the information of all time 
derivatives of the Green function at the shell and then evolving it for later times. 
The time derivatives of the Green function at the shell is given in terms of a recursion 
formula. Using this method we obtain analytic results for short time thermalization 
of the Green function. We show that the late time behaviour of the Green function is 
determined by the hrst quasinormal mode. We then implement the method numer¬ 
ically. As applications of this method we study the thermalization of the retarded 
time dependent Green function corresponding to a minimally coupled scalar in the 
AdS'a and AdS^ thin Vaidya shells. We see that as expected the late time behaviour 
is determined by the hrst quasinormal mode. We apply the method to study the late 
time behaviour of the shear vector mode in AdS^ Vaidya shell. At small momentum 
the corresponding time dependent Green function is expected to relax to equilibrium 
by the shear hydrodynamic mode. Using this we obtain the universal ratio of the 
shear viscosity to entropy density from a time dependent process. 
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1 Introduction 

The question of how a highly excited state in a quantum system relaxes to equilib¬ 
rium or how a quantum system relaxes when one of the parameters describing its 
Hamiltonian is quenched is of phenomenological interest. The relaxation of the state 
produced initially by high energy nuclear collisions from a highly excited state to 
quark-gluon plasma is an important question in the RHIG experiments [1]. Similarly 
the question of quenching of a quantum system can be studied experimentally in cold 
atoms where the coupling of an interacting system can be tuned to almost any value 
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and in short times scales [2-5]. There have been varions approaches, both analytical 
and nnmerical, developed to stndy these qnestions for a variety of qnantnm systems. 
When the qnestion of thermalization is asked in qnantnm field theories which admit 
a holographic dnal, the gange/gravity correspondence links the qnestion of approach 
to eqnilibrinm in the field theory to the formation of a black hole in the bnlk. See 
[6-10] for early stndies pnrsning this idea. 

More recently motivated by snccess of the gange/gravity dnality to describe near 
eqnilibrinm physics and hydrodynamic behavionr in strongly conpled field theories 
this qnestion has received renewed attention which has resnlted in more qnantitative 
nnderstanding [11-17]. In [13] the excited state in the qnantnm held theory was 
created by a translational invariant pertnrbation along the bonndary of a minimally 
conpled scalar held . This pertnrbation lasted for a short dnration of time. It was 
shown by solving the bnlk eqnations that for a small amplitnde of the pertnrbation, 
the metric ontside the in-falling shell of matter is that of a black brane at the leading 
order. This resnlt enconraged snbseqnent anthors to model the collapse to a black 
hole by a homogenons in-falling shell of matter [18-24]. 

For dehniteness we consider the thin shell model of collapse [21] given by the 
following metric 


ds'^ 

e{v) 


4 [-(1 - e{v)z^)dv^ - 2dzdv + dx^] , 

Z2 

f 0, for n < 0, 

\ 1, for n > 0. 


( 1 . 1 ) 


Here z refers to the radial co-ordinate, the bonndary is at ^ = 0. x = x^, ■ ■ ■ x'^~^ 
are the spatial co-ordinates at the bonndary. The metric for n < 0 can be seen to be 
that of AdSd+i nsing the following co-ordinate transformation 


V = t — z. 


( 1 . 2 ) 


While the metric for n > 0 rednces to that of the black brane in AdSd+i nnder the 
co-ordinate transformation 

dv = dt -(1.3) 

1- z^ ^ ^ 

From these co-ordinate transformations, we see that v coincides with time t at the 

bonndary. We have chosen to work with nnits in which the radins of AdS^+i is nnity. 

The radins of the horizon is also nnity. The Penrose diagram of the collapse is given 

in fignre 1. 

In [21] the stndy of how probes snch as two point fnnctions, Wilson loops and the 
entanglement entropy [25] thermalize in the thin shell Vaidya collapse was initiated. 
This stndy was mainly conhned to the saddle point approximation of the probes. 
The behavionr of the probes were characterized in terms of their minimal geometric 
volnme. In [23] the analysis was extended to stndy the two point fnnction of operators 
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singularity 



Figure 1: The Penrose diagram of collapse in AdS Vaidya. 


dual to a minimally coupled massive scalar beyond the geodesic approximation for the 
case of the AdS^ Vaidya shell. The retarded Green function GR{t 2 ,ti'i k) with time 
ti < 0,^2 > 0, before and after the collapse of the shell was evaluated numerically. 
Translational invariance in the spatial direction of the collapsing shell (1.1) enabled 
the characterization of the Green function in the Fourier k space corresponding to 
the spatial directions. The analysis was numerical and it showed that the relaxation 
of the Green function is determined by the hrst quasinormal mode. In [24] the 
study of thermalizing Green functions was extended to fermions in the AdS^ Vaidya 
shell. The analysis was again done numerically. Some analytic properties of the time 
dependent scalar Green functions was studied in [26]. 

In this paper we develop a new method to evaluate the retarded two point 
function GR(t 2 , ti, k) of an operator in the dual theory corresponding to the collapsing 
Vaidya AdS shell in the bulk. The method is general and can be implemented in 
arbitrary dimensions and for arbitrary types of fields in the bulk. The method relies 
on performing the matching of the wave functions corresponding to the dual fields 
before and after the shell term by term in the expansion of the radial co-ordinate z 
of (1.1). We will show that this enables the determination of all the time derivatives 
of Gr just after the collapse, u = O’*" of the shell. It is then possible to evolve the 
Green function to an arbitrary future time t 2 - On implementing this method we see 
that to obtain information of more and more higher derivatives of the Green function 
one needs to perform the matching of the wave functions of the bulk fields to higher 
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powers in the radial coordinate 2 ;. This implies that one needs the information of 
the wave functions of these helds closer to the horizon to obtain long time behaviour 
of the Green function. It also implies that short time behaviour of the two point 
function after the collapse can be determined analytically from the near boundary 
information of the wave functions. Using the fact that the Green function at long 
time is determined by the near horizon behaviour of the wave functions we show 
that the relaxation of the Green function to equilibrium is determined by the first 
quasi-normal mode of the dual bulk held corresponding to the operator of interest 
in the black hole background. 

We implement the method numerically and re-visit the case of the minimally 
coupled massive scalar in the AdS^ Vaidya shell. We reproduce the results of [23]. 
We then study the case of the minimally coupled massless scalar in AdS^ Vaidya 
shell. The Green function corresponding to this scalar is the retarded two point 
function of the spin 2 part of the stress tensor. We show that this Green function 
relaxes by the hrst quasi-normal mode which was determined numerically in [27]. 
We then study the vector fluctuations of the metric and obtain the Green function 
of the vector part of the stress tensor. For small momentum k, it is known [28] that 
this mode admits a hydrodynamic quasi-normal mode which obeys the dispersion 
relation given by 

where rj is the shear viscosity, s the entropy density and T the temperature of the 
fluid. We show the time dependent Green function corresponding to the vector 
fluctuations of the metric relaxes to equilibrium at small k by the hydrodynamic 
mode. Using this we determine the universal ratio of shear viscosity to entropy 
density from a time dependent process. 

This paper is organized as follows. In the next section we detail the new method 
developed in this paper to evaluate the retarded Green function in collapsing AdS 
Vaidya thin shell backgrounds. We see that the method results in a recursion formula 
for the derivatives of the Green function just after the collapse of the shell. In section 
3 we apply the method to obtain the Green function of the operator dual to the 
minimally coupled scalar in AdS^ Vaidya shell. In section 4 we show that the long 
time behaviour of the Green function is determined by the first quasi-normal mode. 
This is first done for the case of AdS^ Vaidya for which wave functions before and 
after the collapse of the shell are known exactly. Then the argument is extended in 
general for any Green function. In section 5 we turn to the case of AdS^ Vaidya. 
We first study the thermalization of the shear correlator of the stress tensor by 
solving the minimally coupled scalar in AdS^ Vaidya. We then examine the vector 
perturbations of the metric to evaluate the two point function of the spin-1 part of the 
stress tensor and show that it relaxes at small momentum by the shear hydrodynamic 
mode. Section 6 contains the conclusions. Appendices A to D deal with technical 
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details required for the analysis in the paper. Appendix E describes the Mathematica 
notebooks which enable the evaluation of the Green functions numerically using the 
recursive method developed in this paper. 

2 Recursion method for time dependent Green functions 

In this section we will outline the general method to obtain the retarded Green 
function in the thin shell Vaidya AdS geometry. The method is general, and can be 
applied to any held in the AdSd+i Vaidya geometry given in (1.1). For dehniteness 
let us focus on the minimally coupled scalar of mass m. The differential equation 
obeyed by 0 is given by 

h{v, Z)dl(j) + ^ 

where 

h{v,z) = l-e{v)z‘^, (2.2) 

and k is the Fourier conjugate of direction 

The solution in the region u < 0 corresponding to before the formation 
of the black hole admits a closed form in terms of Bessel functions. The analytical 
solution will be explicitly discussed in the examples we will consider subsequently. 
The solution has the form 

k, z] ti) = z^*J{v — ti^k^z)^ (2.3) 

where A+ > A_ are the two solutions of the equation 

A(A — d) = m?. (2.4) 

The solution we choose satishes the boundary condition 

/c, 2;; ti) z^-5{v - G) + ..., (2.5) 

with ti < 0. This is necessary to obtain the retarded Green function. In (2.3) note 
that due to time translational symmetry for u < 0, the wave function just depends on 
the combination v — ti. In the black hole region u > 0, the bulk equations of motion 
usually do not admit a closed form solution. But the solution in the frequency u and 
momentum k domain can be constructed in terms of a Frobenius series around the 
boundary z = 0. Then the most general solution in the time domain u > 0 can be 
obtained by taking Fourier transform of the two independent solutions obtained by 
the Frobenius method with respect to the frequency, u. We write this as 

/ OO 

duje~'^‘^'^ ^ [z^+C{u, k)An{uj, k)z^ + z‘^~D{u, k)Bn{uj, k)z"'] . 

n=0 

( 2 . 6 ) 
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Note that here we have assumed that the roots of the indicial equation of (2.1), 
A+, A_ do not differ from each other by an integer. The discussion can be carried 
out for the case when the roots differ by an integer but as we will see that we will only 
need the less singular solution which falls of as z‘^+ to construct the retarded Green 
function. C{uj,k) and D{u,k) in (2.6) are functions which must be determined by 
continuity at n = 0. Note that the coefficients of the differential equation (2.1) are 
discontinuous, but the discontinuity is hnite across n = 0, therefore the solution 0 
must be continuous at n = 0. Thus we have 


lAdS 


{v = 0,k, = (j) {v = 0,k,z). 


(2.7) 


We impose continuity by equating each term of the power series in ^ about the 
boundary. Thus we expand both sides of (2.7) in powers of 2 ; and obtain the equation 


+ ( 2 . 8 ) 

poo ^ 

= / duj [z^+C{oj)An{oj)z'^ + Z^~ D{oj)Bn{oj)z'^) . 


n=0 


n=0 


We have suppressed the dependence of Jn on ti,k and the dependence of C{u), An, 
D{uj),Bn on k to un-clutter the equations. Equating the coefficients of z” in the 
terms proportional to z^+ we obtain 


Jrt. 


du C{uj)An{oj)- 


(2,9) 


From the examples considered in the paper, it is seen that An{oj) is an n-th order 
polynomial in ca Therefore we write An as 


An — ^ ^ 0,n^^ ■ 
j=0 

Substituting this expansion in (2.9) we obtain 

^ poo 

/ dujC{u)u^. 

j=o 


( 2 . 10 ) 


( 2 . 11 ) 


Let us dehne the j-th moment of C (cj) as. 


M,= 


du C{u)u)^. 


( 2 . 12 ) 


^This can be seen using the recursion formula obtained during the construction of the Frobenius 
series solution of (2.1). 
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(2.13) 


Substituting this in equation (2.11), we can rewrite it as, 

n 

■In = 

3=0 


This equation can be inverted to obtain the moments, Mj, which contain information 
about C{u). Note that a{^ are known from Frobenius series solution cj)^^. Therefore 
we obtain 


o-iMjj ; n > 0, (2.14) 

Mo = 4- (2.15) 

% 

Knowledge of all the moments Mj is sufficient to construct the retarded Green 
function. To see this consider the near boundary behaviour of the field 0. From 
equations (2.5), (2.6) and (2.8), the near boundary behaviour of the held for u > 0 
is given by 


Mn — 

aH: 


J n 


n-1 

E 

1=0 


lBH 


{v,tuk)= I duje-^^^C{uj)z^Ao + --- + 0{z^-) 


(2.16) 


Together with the boundary condition (2.5), the AdS/CFT recipe for the retarded 
Green function [29] states that it is given by 

/ OO 

due-^‘^'^C{u). (2.17) 

•OO 

Here we have ignored overall proportionality constants in the Green function to 
simplify the discussion. Expanding the exponential as Taylor series, and using the 
dehnition of the moments of C{u) we obtain 


Gr{v) = 

n=0 




(2.18) 


Thus the knowledge of all the moments of C{u) is sufficient to construct the Green 
function. Note that knowledge of D{u) in (2.6) is not necessary. We call this the 
recursion method to obtain the Green function since each term in (2.18) is given 
recursively from the knowledge of the lower moments using (2.14). It allows for the 
construction of the the Green function as a power series in time for u > 0. The 
information of ti is present in and all terms depend on momentum k. There is 
another way to view this construction of the retarded Green function. Note that 
the j-th moment of C{u) is essentially the j-th derivative of the Green function, 
evaluated at u = O’*", 


dv^ 


11=0+ 


du{—iujyC{u) = {—iyMj. 


(2.19) 
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Thus this method determines the Green function for n > 0 from the knowledge of all 
its time derivatives at n = O'*". 

We now make some general properties of this method of determining the retarded 
Green function. Since the method relies on construction of using the Frobenius 
series we can apply it in general to all black hole backgrounds even if closed form 
solutions do not exist. The method can be applied even if the background is only 
known in terms of a power series in about the boundary for example in hnite 
temperature versions of RG flow solutions [30] . Note that for short time development 
of the Green function for n > 0, knowledge of only a few moments is needed. This 
implies from (2.18) and (2.14) we need the knowledge of wave functions before and 
after n = 0 close to the boundary. That is we need the knowledge of the wave 
functions to a few powers of 2 ;. Thus short time development of the Green function 
can be written down analytically by obtaining a few moments. However for the long 
time behaviour of the Green function we need to know a large number of moments. 
Again from (2.14), this implies we need the information of the wave function for large 
powers of ^ which in turn implies that we need the behaviour close to the horizon. 
This hts with the general intuition that long time behaviour is controlled by the 
behaviour near the horizon. It is this property which will enable us to prove that the 
long time behaviour of the Green function is determined by the quasinormal mode in 
section 4. Finally, it will turn out that even moments are real and odd moments are 
purely imaginary. This is because the coefficients for odd j are imaginary. This 
is easily seen due to the fact that each power of u comes with a factor of i. Then 
from (2.14) it is easy to see that even moments are real and odd moments are purely 
imaginary. This then ensures that the Green function given by (2.18) is real. 


3 Thermalization in AdS^, Vaidya 

In this section we will implement the method developed in section 2 for the case 
of thin shell Vaidya metric in AdS^. We study the thermalization of the retarded 
Green function of the operator corresponding to the minimally coupled scalar of mass 
m. We will show that using our method we reproduce the results of [23]. We also 
determine a few low moments analytically to obtain the short time behaviour of the 
Green function. A simplihcation that occurs for the case of the AdS^ is that the 
solutions of the minimally coupled scalar in the BTZ black hole are known in closed 
form in terms of hypergeometric functions [31]. These solutions will provide the 
initial starting point in our argument to demonstrate that the long time behaviour 
of the Green function is determined by the lowest quasi-normal mode. 


The thin-shell AdS-Vaidya metric, in 3 dimensions is given by 



h{z, v) 


h{z,v)dv^ — 
= 1 - e{v)z‘^. 


2dvdz + dcp^], 


(3.1) 


The co-ordinate 0 parametrizes the spatial direction of the held theory and we assume 
that it is not compact. The scalar held equation in this metric is given by 

- 2dyd^^^ + (^d^h - 

3.1 Scalar wave functions in AdS^ Vaidya shell 

We solve the minimally coupled massive scalar equation given in (3.2). We discuss 
the solutions before and after the collapse of the shell below. 



+ (3.2) 


Solution for n < 0 


The solution for the scalar held for n < 0 is given by [23] 


$ 


AdS/ 


= c 


ti,k,z) 

6{v 


ti)z^ 


[(n - fi)2 2{v - ti)z] 


\k\ 


u+ 




_ 2V>'0F 

r(0 


(3.3) 

(l^l\/(^ - fi)2 + 2(v - ti)z^ , 


(3.4) 


Note that this solution is written down in mixed Fourier space {t, k). Here A± are 
solutions to the quadratic equation A(A — 2) = and are given by 


A± = 1 ± z/, u = y/l + m?- 


(3.5) 


This solution satishes the required boundary condition discussed in (2.5) which is 
required to evaluate the retarded Green function, 

— ti, k,z) = 6{v — ti)z^~ + _ (3.6) 


To show this we hrst take the the limit ^ 0 in (3.3) keeping v — ti 0. We see 

that there is no term proportional to z^-. Now we take v —)■ H, the wave function 
then reduces to 




0{v — tl)2^/n z^+ 

r(-// + l)r(i') [(tf - ij)2 + 2(u 



(3.7) 


This solution certainly diverges in the limit n —)■ ti. All what one needs to show is 
that the expression in (3.7) is a representation of the delta function times z^~. To 
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demonstrate this we perform the integral over v as follows 




2^ 


A_ 


= ^ 


r(-^ + i)rM 

A_ 


dv 

— 9{v - ti) 


[(^ + 1)2 - 1]-+^ 


(3.8) 


Therefore we conclude that the solution in (3.3) satishes the required boundary 
condition given in (3.6). 

Solution for u > 0 


The strategy to obtain a closed form solution for u > 0 is as follows. Note that 
under the transformation 

‘ ^(^) ■ 

the metric given in (3.1) for u > 0 reduces to that of the BTZ black hole in Poincare 
coordinates. This is given by 


^ (-(1 - z^)de + I. 


(3.10) 


Now the solutions to the minimally coupled scalar in the BTZ black hole is known 
[31]. The two independent solutions are given by 


= e 






[l-zT^z 


(3.11) 


- ^))> + k)), z^ ] , 






2n ^A_ 


(1 — Z^) 2 Z 


xF (-{A_ + i{u + k)), -(A_ + i{u — k)), 1 + iu, 1 — z^ 

y 2 2 

Note that these two independent solutions reduce to the ingoing and outgoing Fourier 
modes at the horizon. We can now obtain the solution in the coordinates {z, v, 0) by 
performing the substitution given in (3.9). This leads to the following independent 
solutions for the Vaidya metric in the region u > 0, 


= e 


_ —iu}v-\-ik4> 


(1 + zY^z 


iu) ^A_ 


(3.12) 


( ^(A- - i{u - fc)), ^(A_ - + k)), ] , 




_ —iLOv-\-ik(f) 


(1 - zy‘^z 


A_ 


' 1 1 

xF ( + i{uj + k)), 2 ^^- + ~ I + iuj,l - 


- 10 - 







Now the above solutions admit an expansion around the horizon z = 1, however we 
have seen in section 2, to obtain the Green function we need an expansion around 
the bonndary z = 0. This can be achieved by using the transformation properties 
of the hypergeometric fnnctions. We have performed the reqnired transformation in 
appendix A. This results in the following independent solntions, 

(1 + zy^ X (3.13) 

^ ~ ~ ~ ’ 

F Q(^+ - + F)), ^(A+ - i{u - k)), A+, z^'^ . 


It can be explicitly verified that the above solntions satisfy the eqnations of 
motion given in (3.2) for n > 0. The partial Fonrier transform of this solntion with 
respect to cj is given by 


+ 


poo 

z^- / due-^‘^'^Ci{u){l + z)-^^F 


’ —OO 

roo 


ckje-‘“"C 2 {ui){l + z)-‘“F 


(3.14) 

Q(A- - i{uj - k)), ^(A_ - i{uj + k)), A_, 
Q(A+ - i{,uj + k)), ^(A+ - i{LJ - k)), A+, z‘^^ . 


The above solution admits an expansion around the bonndary z = 0. Therefore we 
have written the solution in the black hole region in the reqnired form given in (2.6). 
Note that in the case of the BTZ black hole, the Frobenius expansion aronnd the 
bonndary can be written in closed form. 


3.2 Matching at n = 0 and Green function 

We follow the general procednre discussed in section 2 to constrnct the time depen¬ 
dent Green fnnction. To do this we hrst obtain the moments of the fnnction ^ 2 ( 0 ;) 
by matching the wave fnnction k,z) in (3.3) and ^^'^^{v,k,z) in (3.13) at 

n = 0. The moments of ^ 2 ( 0 ;) are determined by comparing powers of x in the terms 
proportional to z^+. We will demonstrate this procedure explicitly and obtain mo¬ 
ments np to the second order. We expand the LHS of (3.3) and the coefficient 
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of (3.13) to quadratic order in 2 ; to obtain the following equation 


ce{-h) 












dz 


\z=0 


+ 


~2 dz^ 


/ dujC 2 {oo){l — iooz + iu{iu + 1) — + ...)x 

-CX) 

- Ul + V - i{uj + k))Ul + V - i{uj - k)) „ . 

X (1 + ^^+ ...) 


duC 2 {u)){l — iuz + 


1 + z/ 

1 + u u? 

y 


\ 2 \ 

z^+ ...). 


4(1 + z/) 


(3.15) 


Comparing the terms we can read out the moments to quadratic order. These are 
given by 


Mo 

Ml 

M 2 


' —00 
roo 


duC2(uj) = A^J_^_i(|A;||ti|), 

dujwC2{uj) = —iNkJ_^_3{\k\\ti\), 


(3.16) 


' —00 

POO 


du w‘^C2{oo) = N 


(1 + z/)^ 2 2(1 + z/)(l + 2z/) 


3 + 2z/ 
2k{l + z/ 
U 


+ r- 


''1 




j 


-v+ 




here, N = C6*(—fi) As discussed in section 2, see equation (2.18), the 

retarded Green function is given by ^ 


G{v) = Mo + ivMi - ^^ v^M 2 H-. (3-17) 

Thus it is clear that the short time expansions of the Green function can be easily 
obtained. 

It is interesting to take the fc = 0 limit of the moments and construct the short 
time expansion of the Green function. In this limit the Bessel function reduces to a 

^ Since we are interested only in the time dependence we are ignoring overall proportionality 
constants in the Green function. 
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rational function and the first few moments are given by 


Mo 

Ml 

M2 


C 2^^+! 


ff+ir(-z/+ D’ 

C 2^^+! -i2{-u+l) 


C 2'^+! 1 +1/ 

ff+ir(-z/+ i) (3 + 2z/) 


1 + z/ — To 2z/)(3 + 2z/) 


''1 


(3.18) 


Note that it is clear that the Green function does not vanish identically for fc = 0 for 
arbitrary t 2 ,v, since these moments do not vanish. 


3.3 Recursive numerical construction of the Green function 

It is easy to set up an algorithm in Mathematica to evaluate the moments recursively 
as discussed in section 2. This algorithm is used to evaluate 56 moments of the 
function ^ 2 ( 0 ;). From (2.18) we can construct the Green function to 0{v^^) ^ To 
improve accuracy we then approximate the Green function using the (28|28)th Fade 
approximant. The results for the Green function are given in the 3 figures which we 
will discuss. 

Figure 2, shows the thermalizing Green function (solid blue curve), the vacuum 
Green function (dot-dashed red curve), and the thermal Green function (green dashed 
curve), as a function of future time v. The thermalizing Green function starts close 
to the vacuum Green function, for small time, however it deviates away from the 
vacuum Green function within one horizon time. At large time, the thermalizing 
Green function approaches the thermal one. At n = 0, the thermal Green function 
starts at a different value than the thermalizing and the vacuum Green functions. 
The thermal Green function is plotted from (B.2) derived in the appendix B. Figures 
2 and 3, are plotted for specific values, z/ = |, = | and ti = —1.7. 

Figure 3 is the logarithmic plot of the absolute value of the thermalizing Green 
function (solid blue curve), and the imaginary part of the lowest quasinormal mode 
(dot-dashed red curve), It is seen that for large time, the decay of the ther¬ 

malizing Green function is given by this lowest quasinormal mode. Here, large time 
means time of the order of a few horizon radius, as can be seen from the plot. It 
is important to note that thermalization, i.e. decay of thermalizing Green function 
follows the lowest quasinormal mode, is achieved within a few (~ G(l)) horizon 
radius. The dips in the plot are the points where the thermalizing Green function 
crosses the time axis in figure 2, and indicates the oscillations of this Green function. 
Figures 2 and 3 reproduce those found in [21] by the direct numerical integration of 
the differential equation (3.2). 

^We stopped at this order in moments since we found that for times n = 5 in horizon units the 
Green function evaluated converged to high degree of accuracy. 
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G«(v,?i,k) 



Figure 2: The thermalizing (solid blue curve), vacuum (dot-dashed red curve) and 
thermal (dashed green curves) Green functions are plotted as a function of future 
time, u, for hxed values, ^ = I, k = ^ and ti = —1.7. 

|GK(v,?i,k)| 



Figure 3: The logarithmic plot of the absolute value of the thermalizing Green 
function is given by the solid blue curve, for fixed values, i/ = |, fc = | and ti = —1.7. 
The dot-dashed red line gives the lowest quasinormal mode of the thermal Green 
function. 

Lastly, in figure 4, the thermalizing (red dots) and the thermal (solid blue curve) 
Green function are plotted as function of k. The value of the thermalizing Green 
function close to /c = 0 is non zero, as expected from the earlier discussion in around 
(3.18). Also, the values of the thermalizing and the thermal Green functions are close 
to each other near k = 0. This figure is plotted for specific values, z/ = |,ti = —2 
and t 2 = 5. 

4 Long time behaviour of the Green function 

It is clear from the expression for the Green function given in (2.18), that for obtaining 
the long time behaviour of the Green function we need to evaluate moments for 
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Figure 4; The thermalizing and thermal Green functions as functions of k are given 
by the red dots and solid blue curve, respectively, for fixed values, i/ = |,fi = —2 
and t 2 = 5. From the plot it is seen that values of the thermalizing and thermal 
Green function are close to each other near fc = 0. 

large values of n. Form the recursive relations for the moments in (2.14) we see that 
this can be done if we implement the matching of the wave functions to order 
This in turn implies that we need the knowledge of the less dominant wave function 
in (2.6) dehned by 

OO 

$■•■( 0 ;, k, z) = ^ z^+An{uj, k, z)z"', (4.1) 

n=0 

closer to the horizon. This is because the coefficients of z^ for large values of n will 
be determined by the singular behaviour at the horizon. The behaviour of this wave 
function near the horizon can be determined by the general properties of solutions 
of wave equations in AdS black holes. Using this information we will show that the 
long time behaviour of the Green function is determined by the first quasinormal 
mode. We will first demonstrate this for the case of the minimally coupled scalar in 
AdSs Vaidya. This is easy to do explicitly since, the wave function is known 

in closed form. We will then show in general that the long time behaviour of the 
Green function is determined by the first quasinormal mode of the wave functions in 
the corresponding black hole background. 

4.1 Green function in AdS^ Vaidya 

From the preceding discussion we need to examine the behaviour of the function 
$■*■( 0 ;, k, z) close to the horizon. For the case of the minimally coupled scalar in the 
BTZ black hole this function is known in closed form (3.14) and is given by 

^STZ(+)(^,^,^) ^ + k)),^{A+ - i{u - k)), A+, z^^ . 

(4.2) 
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The near horizon limit of this wave function can be easily determined from the 
properties of the hypergeometric function given in (A.l). We see that the singular 
behaviour near z = 1 is given by 




r(A+)r(-zn;) 


r|(A+ - i{oj - fc))r(|(A+ - + k)) 


(4.3) 


This equation will serve as the starting point of obtaining the long term behaviour 
of the Green function. It is important to note that the mode 'h’*' vanishes near the 
horizon when the frequency is given by 


ujY = —iA_|_ — 2m ± A;, n = 0,1, 2, • • • . 


(4.4) 


These are the quasinormal modes of the minimally coupled scalar in the BTZ black 
hole [32]. Also note that these modes are located in the lower half cn-plane. We now 
impose the matching condition at n = 0. This will lead to the following equation 


C 


\k\ 


u+i 


-J_ 


/ n , , N (2v + l) 

{tl + 2\ti\z) 4 

'OO 

duC 2 {uj){l — z)~ 


tf + 2|fi 1^) 


(4.5) 


r(A+)r(-za;) 


r|(A+ - i(uj - fc))r(|(A+ - i(uj + k)) 


We expect this approximation near the horizon to estimate the behaviour of the 
moments suitable to obtain long time behaviour of the Green function. We can now 
solve for C 2 (cz) by substituting 


J/ = ltl(l -z), 


(4.6) 


and multiplying both sides of the equation in (4.5) by e and formally integrating 
over y from —oo to oo This leads to the following expression for 6 * 2 (ca) 


^2 


.cz 


/(ca, k) 


ri(A+ - i(u - k))r(l(A+ - z(u + k)) 

r(A+)r(-^a;) 


(4.7) 


where 


I(u, k) = — 


C 


dy 


\k\ 


^+5 


{tf + 2\ti\z)^"4^ 




(4.8) 


Finally the Green function is given by 

/ OO 

due-^^^C2{oo) = 

■OO 

.-i..„rl(A+ - i(u) - «:))r(l(A+ - i(u + k)) 


(4.9) 


due 




k). 


^Though the range of y is restricted from —oo to 0, we are extending this range formally to 
obtain C 2 - 
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Note that the factor multiplying I{u,k) in (4.9) is a function of u such that for 
to —)■ —ioo it behaves as 


lim H{u, k) 

U)^ — iOQ 


= lim 


r|(A+ - i{uj - fc))r(4(A+ - i{uj + k)) 


r(A+)r(-za;) 




h'H- ? 


(4.10) 


/(cu, k) is essentially a Fourier transform of the Bessel function. We assume that this 
is analytic in the lower half cu-plane and grows at the most exponentially in to as 
00 —)■ —ioo, given by 

lim i{oj,k) (4.11) 

LO^ — iOQ 

Here M > 0 is a fixed constant. With these assumptions and the behaviour in (4.10) 
we see that integrand in (4.9) goes to zero in the limit oo —>■ —ioo for sufficiently 
large but fixed u > 0. This is because of the exponentially dying term in 

the integrand of (4.9). Therefore the integral can be performed by completing the 
contour in the lower half cu-plane. Then the integral localizes to a sum over the poles 
of the Gamma functions. The poles of the Gamma functions are at the quasinormal 
modes given by (4.4). The result of the integral then reduces to 

OO 

G{k,v,ti) r-u ^ fc)(—27ri Residue^=t^a if(a;, fc)). (4.12) 

n=0,a=it 

It is now clear from (4.12) that the long time behaviour of the Green function is 
determined by the lowest quasinormal mode as we have seen in our explicit numerical 
evaluation of the Green function. The decay of the Green function is controlled by 
the imaginary part of the lowest quasinormal mode and the period of oscillations is 
determined by the real part of the lowest quasinormal mode. Note also from (4.12) 
we see that since the expression involves a sum over all the quasinormal modes, the 
rough time scale over which the lowest quasinormal modes takes over is of the order 
of a few horizon times. This is also clearly seen in the numerical evaluation of the 
Green function. It is important to note that the starting point of the analysis was the 
behaviour of near the horizon given in (4.3), which vanished at frequencies 

determined by the quasinormal mode spectrum in the lower half cu-plane. 


4.2 Green function in AdSd+i Vaidya 

Using the intuition gained by the explicit solutions of the minimally coupled scalar 
in the BTZ black hole we now generalize the discussion. We show that the long time 
behaviour of the retarded Green function in arbitrary AdS^+i Vaidya background is 
determined by the lowest quasinormal mode of the corresponding bulk field. 

Gonsider the differential equation given in (2.1) for the minimally coupled scalar 
in the AdSd+i Vaidya for u > 0. Substituting 

0(u, k, z) = ^uj{k, (4.13) 
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we obtain 


^ - 2d^^uj ) - 


dz 


d-l 




m 


+ j $^ = 0 . 

(4.14) 

This equation has two regular singular points, z = 0 and z = 1, corresponding to 
the boundary and the horizon of the black hole. One can set up a Frobenius solution 
around either z = 0 or z = 1. Let the two independent solutions around = 0 be 
and <F“. From the indicial equation for the expansion at z = 0 we know that 
these solutions behave as follows ^ 


lim^^fz, k) ^ z‘^+, lim^ (z, k) z‘^~. (4-15) 

z—^O z—^0 

Therefore <F“ is the dominant singular mode at the boundary. Similarly one can set 
up an expansion around z = 1, let the two independent solutions around z = 1 be 
It is easy to see from the indicial equation around z = 1, that their 
behaviour near 2 ; = 1 is given by 

lim ~ lim K"(z,k) ~ (1 - z)°. (4.16) 

2->l 2->l 

The reason we have labeled these solutions as in and out is because they correspond 
to the ingoing and outgoing solutions when these wave functions are transformed to 
the t, z coordinates. To see this note that the coordinate transformation near z = 1 
can be obtained by integrating (1.3). This is given by 

n ~ t + - ln(l — z). (4-17) 

(Jj 

Substituting this coordinate transformation we obtain the corresponding wave func¬ 
tions in the t, z coordinates 

rf{v,z) = e-^‘^^¥j:{z,k) ~ e-*‘^*(l - ^)^, (4.18) 

<Pl^iv,z) = ~ e-‘"*(l - z)-^. 

Now we have 2 sets of linearly independent solutions and {*F^, 4)°“*}. 

Therefore by the uniqueness theorem of second order ordinary linear differential equa¬ 
tions we should be able to express one set in terms of the other as linear combinations. 
Let us write 


$r(A:,^) = Mn{u,k)<!>+{k,z) + Mu{u,k)<!>Z{k,z), (4.19) 

$L"(fc,;2) = M2i{u,k)^^{k,z) + M22{u,k)^Zik,z). 

®We have assumed that the roots of the indicial equation do not differ by an integer. The 
discussion can be easily generalized in case they do. 
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Let us now use the definition of quasinormal modes to obtain some information of the 
coefficients M 12 and M 22 - We have seen that correspond to the ingoing and 

outgoing modes in the f, ^ coordinates. Now by definition quasinormal modes are 
those values of frequencies for which these modes obey Dirichlet boundary conditions 
at the horizon. This implies that the more dominant mode proportional to should 
vanish at these frequencies. Thus we have the equation 


Mi2(a;r,^) = 0, M22(a;r,^) = 0, 


(4.20) 


where 0 ;°“* are the quasinormal frequencies which he in the upper half (u-plane for 
the outgoing modes and are the quasinormal frequencies which he in the lower 
half ca-plane for the ingoing modes [32, 33]. Let us invert the equation (4.19), we 
obtain 


\ ^ ( M22 -Mi 2 \ 

<h-(fc,z )) detM V-M 21 Mn ) V ) ' 


Note that the inverse exists that is detM ^ 0. This is because and <h,, can be 
written as linear combinations of and We can now easily read out the 

near horizon behaviour of the solution from the above equation. The singular 
behaviour of is given by 


- MMoJ,k)t‘:{k,z )), (4.22) 


’S>:{k,z)\ 


z-s-1 


detM 


M22{uj, k){l - zy d . 


Here we have used the near horizon behaviour given in (4.16). Now M 22 {oj,k) van¬ 
ishes at UJ = with zeros in the lower half plane. This is the general form of 
the equation (4.3) seen explicitly for the case of the BTZ black hole for black holes 
in AdSd+i- Note that we arrived at this result from the general definition of quasi¬ 
normal modes. From this point onwards we can follow the rest of the argument in 
section 4.1 to arrive at the conclusion that the long time behaviour is determined 
by the leading quasinormal mode. In general the gaps between quasinormal modes 
are of the order of horizon scales. Therefore we expect the leading behaviour to set 
in order of a few horizon times. Though here we have used the minimally coupled 
scalar to demonstrate our argument for simplicity, the analysis can be carried out 
for other fields. The steps involved will result in similar equations. We will see this 
explicitly for the vector fluctuations of the metric in the subsequent section. 

This concludes our general argument of why the long time behaviour is deter¬ 
mined by the leading quasinormal mode. Thus any retarded correlator which is used 
to probe the onset of thermalization caused due to injection of energy at an instance 
of time in the field theory will decay with the time scale set by the first quasinormal 
mode. In [34] it was conjectured that thermalization time scales in the field the¬ 
ory are determined by quasinormal modes. We have shown that the time dependent 
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Green functions considered in this paper provides an explicit realization of this state¬ 
ment. From our general argument we expect this to be true for correlators whose 
lowest quasinormal modes are determined by hydrodynamics if the probing momen¬ 
tum scales are sufficiently small. In the next section we will verify this expectation 
for stress tensor correlators in A/" = 4 Yang-Mills whose lowest quasinormal mode is 
sensitive to the shear viscosity. 

5 Thermalization in AdS^ Vaidya 

In this section we study the thermalization of the retarded two point functions in 
the thin shell AdS^ Vaidya geometry. We hrst consider the equation satisfied by the 
spin 2 metric ffuctuation This ffuctuation, perpendicular to the momentum k 

which is along the direction, forms the shear spin 2 mode. Evaluating the retarded 
Green function holographically for this mode provides information of the retarded 
two point function of the stress tensor , for the strongly coupled A/" = 4 

Yang-Mills The lowest quasinormal mode of shear fluctuations of the metric has 
been studied earlier in [27] numerically. We see our results are consistent with this 
earlier calculation. The more interesting metric fluctuation to consider is It is 

known that this mode admits a hydrodynamic quasinormal mode in the AdS^ black 
hole for small momentum [28] . This mode appears as a pole in the thermal two point 
function of the stress tensor . We will implement the recursion method 

numerically and evaluate the thermalizing Green function for various values of small 
momentum k. From the long time behaviour of the Green function we show that this 
Green function relaxes by the hydrodynamic quasinormal mode. This enables us to 
read out the universal shear viscosity to entropy density ratio from a time dependent 
process. 

The thin shell AdF-Vaidya metric in 5 dimensions is given by, 

ds"^ = \ \—h(v, z)dv^ — 2dzdv + dT^] , (5.1) 

h{y,z) = 1 — 9{y)z'^, (5.2) 

here, T is a 3 vector, with components, x®, where i = 1,2, 3. The a:®’s parametrize 
the spatial directions in the dual field theory. The transformation which reduces the 
n > 0 part of the above metric to the AdS^ black hole metric in Poincare co-ordinates 
is given by 

d z 

dv = dt — - -(5.3) 

-L Z 

^Metric perturbations for spherical shell collapse was studied in [35]. Non-equilibrium Green 
function corresponding to the shear correlator was studied earlier by a complementary approach by 
examining backgrounds perturbed by hydrodynamic modes [36, 37] 
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Figure 5: Figure (a) shows the thermalizing (solid blue curve), and vacuum (dot- 
dashed red curve) Green functions for the shear mode. Figure (b) shows the loga¬ 
rithmic plot of the absolute value of the thermalizing Green function. Plots are for 
values = —15 and k = 10. 


In these coordinates we obtain the planar black hole metric in AdS^, for n > 0, 


ds^ 


1 


— (1 — z^)dt^ -I- 


dz^ 

1-z^ 


- 1 - dx^ 


(5.4) 


Here again the radius of AdS as well as the black hole has been set to unity. The 
temperature of the black hole is then given by 

T=1 (5.5) 

71 

To study two point functions of the stress energy tensor of the boundary field 
theory we consider small perturbations to the background metric, , of (5.1), 
dfiu = gjS + The differential equations satished by the metric perturbations 
are obtained by linearizing the Einstein’s equations which is given by 


'klfiu =-4:g^u- (5.6) 

Here is the Ricci tensor and the value of cosmological constant. A, has been set 
to A = —6, in units of AdS radius. The linearized Einstein’s equations are given by 

K‘t = -4V. (5.7) 

where is the linearized Ricci curvature. 


5.1 Spin 2 metric perturbations 

First, we consider the case when only the scalar mode of metric perturbation, 
is turned on. We substitute the ansatz ^ z) = in (5.7). 
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Here, the momentum is chosen to be along the direction. After the redehnition 
$( 2 ;) = z‘^h^ix 2 {z) we obtain the following equation from (5.7), 

u < 0 : z^^" - z(3 - 2iu;z)^' - (k^z^ + Sicjz)^ = 0, (5.8) 

u > 0 : z^(l — z^)^" — z(3 + z^ — 2iujz)^' — (k^z^ + 3iu;z)^ = 0, (5.9) 

where the prime means derivative with respect to z. Note that as expected these are 
the equations of the minimally coupled massless scalar in the AdS^ background for 
u < 0 and the AdS^ black hole background for u > 0. 


Solution for u < 0 


In appendix C we have obtained the solution which satisfies the boundary condi¬ 
tion (2.5) for a minimally coupled massive scalar in AdS^. On setting the parameter 
= 3 or 1 / = 2 in (C.l), the equation reduces to (5.8). Therefore we can read out 
the solution in mixed Fourier space which satisfies the required boundary condition 
from (C.18). This is given by 


k, z] ti) = k, z] ti) = 9{v — ti)2 2 x 

k 


(5.10) 


\/{v - tiY 2{v - ti)z 


Z^J_5 


V - fi)2 2(v - fi)^j 


Solution for u > 0 


Unlike the case of the BTZ black hole, closed form solutions to minimally coupled 
scalar in the AdS^ black holes are not available. However the recursive method we 
developed in section 2 just relies on information of the Frobenius expansion of the 
solution near the boundary. The roots of the indicial equation of (5.9) are given by 
A± = 4, 0. As discussed earlier to obtain the retarded Green function we need to 
obtain the Frobenius solution with the less dominant mode, that is with A+ = 4. 
The recursion relations for the coefficients a^s defined in (2.10), necessary for the 
construction of the solution in the region u > 0, are given by the following equations 
with Oq = 1, 


a 

a 


0 

1 


0 

2 


a 


0 

3 


0, a\ 

k'^ttQ 


= -ial, 
4 = 0) 



0 , 


— i9a2 
21 


—i7a\ 

12 ’ 

4 = 0) 


—i9a\ 

21 


(5.11) 
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and, for n > 4, 


a„ = 


ai = 


ai = 


ai = 


fc^a°2 + n2a°. 


n—4 


n{n + 4) 

-i{2n + 3)a{l\ + + n^al^-A 

n{n + 4) 

-i{2n + 3)a{l\ + 


(5.12) 


; j = l,2,...,n-4, 


n{n + 4) 
-i{2n + 3)a{l\ 


; j = n-3,n-2, 


n{n + 4) 


; j = n-l,n. 


We can now extract the moments of the Green fnnction by snbstitnting the val- 
nes of determined from (5.11),(5.12) into (2.14). The valnes of the coefficients 
Jn in (2.14) can be obtained from the recnrsion relation (D.4) with v = 2. Finally 
the retarded Green fnnction is constrncted nsing (2.18). This procednre is clearly 
algorithmic and can be easily implemented nnmerically. Using Mathematica we im¬ 
plemented this procednre of obtaining the Green fnnction. We have evalnated np to 
200 moments ^ 2 ( 0 ;) and the Green fnnction is constrncted to 0(n^°°) in fntnre time, 
V. We then approximated the Green fnnction by the (1001100) Fade Approximant for 
better accnracy. For the valnes of k stndied, 200 moments were snfficient to obtain 
convergent resnlts. 

Fignre 5a shows the thermalizing Green fnnction (solid bine curve) for fi = —15 
and a particular value of A; = 10. As expected, the thermalizing Green function 
starts close to the vacuum Green function (red dot-dashed curve) at small time, and 
deviates from it within one horizon time. From the logarithmic plot, hgure 5b the 
slope as shown by the dot-dashed red line of the thermalizing Green function (solid 
blue curve) is measured at large time, i.e. at time of the order of a few horizon 
radius. This slope gives the negative imaginary part of the lowest quasinormal mode 
(—Im ca), as shown in section 4. From the gaps between consecutive zeroes of the 
thermalizing Green function, real part of u can be calculated using, Reca = vr/gap. 
The values of real and imaginary part of cn, for several values of k, are listed in 
table 1. These values agree with hgures 5,6 of [27] where the numerical values of the 
real and imaginary part of the quasinormal modes for the minimally coupled scalar 
with A = 4 in AdS^ was obtained for various values of k. To compare our results 
to that of hgures 5,6 of [27], note that we need to perform the following scalings 
Incurs _ 2 qtheirs^ ^ours _ 2 ^theirs^ simple check uote that = 10 corresponds 

to = 5^ looking at their hgure 6 we note that = .8 which corresponds 

to Ima;““'’* = 1.6 which agrees with our result in table 1. Similarly note that for 
small values of the value of is above the 45° line. This is also the 

case from our results in table 1. 
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k 

Re u 

—Im u 

7.0 

8.34 

1.90 

7.5 

8.84 

1.77 

8.0 

9.29 

1.76 

8.5 

9.80 

1.78 

9.0 

10.28 

1.76 

9.5 

10.72 

1.63 

10.0 

11.21 

1.63 

10.5 

11.70 

1.52 

11.0 

12.22 

1.55 


Table 1: Values of the real and imaginary part of the lowest quasinormal modes, 
—Im (u, calculated from the large time behaviour of the thermalizing Green function, 
for several values of momentum, k, and hxed value of ti = —15. 


5.2 Vector metric perturbations and shear viscosity 

It is known that vector metric perturbations in the AdS^ black hole with 

momentum is along the direction admit a hydrodynamic mode at small momentum 
in addition to the usual gapped quasinormal frequencies [28]. The hydrodynamic 
quasinormal mode corresponds to the hydrodynamic pole in the thermal correlator 
This was used to read out the ratio of shear viscosity to entropy density in 
[28]. Using the methods developed in this paper we can evaluate the time dependent 
thermalizing retarded two point function in the AdS^ thin shell Vaidya 

background. From the general analysis of section 4 we expect that the time dependent 
Green function in the AdS^ thin shell Vaidya background should relax to equilibrium 
by the hydrodynamic quasi-normal mode. Therefore from the decay it should be 
possible to read out the ratio of shear viscosity to entropy density from a dynamical 
Green function. In this section we perform this analysis using the methods developed 
in this paper and obtain the universal ratio of shear viscosity to entropy density^. 

We hrst turn on the following metric fluctuations with momentum along the 
direction, 

K^^{t,x^z) = K.,3it,x^z) = (5.13) 

We can obtain the linearized equations of motion from (5.7). Redehning the helds 
as, 

Hy{z) = z^K^i{z), H‘i{z) = z‘^h^i^3{z), (5.14) 

^As far as the authors are aware this is the hrst instance where the universal ratio of shear 
viscosity to entropy density is obtained from a time dependent process. 
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we obtain the following coupled equations in the AdS^ background before the collapse 
of the thin shell, n < 0, 


H” - -H', + ikH'^ = 0, 


(5.15) 


H”+(-- + iuj iJ' - - kuHs = 0, 

HI + ( -- + if' - —fig + ikH'^ - —H^ = 0. 
' z z z 


After the formation of the AdS^ black hole for n > 0, the equations are given by 


if" - -H', + ikH's = 0, 


(5.16) 




^ J zf f zf 


where / = 1 — z^. The above differential equations can be decoupled to give the 
following equations by considering ii' = Pv, 


V < 0 ■. z^p'l + z{—3 + 2iu)z)pl — (—3 + k‘^z‘^ + 3iuz)p^ = 0, (5-17) 

n > 0 : z^{l — z^)pI + z{—3 — + 2iujz)pl + (3 — k‘^z‘^ + 9z^ — 3iojz)p^ = 0. 

(5.18) 

Solution in AdS^: n < 0 

We see that on substituting w? = 0, i.e. i/ = 1 in (C.l), the equation reduces to 
(5.17). Therefore we can read out the solution in mixed Fourier space which satisfies 
the boundary condition (2.5) from (C.18). This is given by 

k, z] h) = k, z] ti) = (5.19) 

6 {v — ti)2~2y/n 

Solution for AdS^ black hole: n > 0 


k 


V — fi)2 + 2{v — ti)z 


z^J-3 (^k^/{v — fi)2 + 2(n — ti)z'^ 


We again use the Frobenius expansion around the boundary z = 0 to solve 
(5.18). We obtain the following recursion relation for for the the root A+ = 3 of 
the indicial equation where Oq = 1, 


a 

a 


0 

1 


0 

2 


a 


0 

3 


0 , 

k'^ttQ 


a^ = 


-la, 


0 


0 , 




0) 

Ci‘2 — 0 , ^‘2 

k‘^a\ — i7a2 


On = 


—iba\ 


15 


= 0 , 


—i7a\ 

15 


(5.20) 
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(a) (b) 


Figure 6: The solid blue curve is the thermalizing Green function. In (a), the 
dot-dashed curve is the vacuum Green function. Figure (b) has the logarithmic plot 
of the absolute value of the Green function, plotted with its late time slope (red 
dot-dashed curve). Plots are for values, = —5 and k = 1.1. 


and, for u > 4, 


k‘^cbn -2 + — 2n — 8 )a °_4 

n{n + 2 ) ’ 

—i{2n + l)a^T\ — 2n — 8 )a {_4 

n{n + 2 ) 

. -i(2n + lXl\ + fcX-2 

ai = -^-; j = n — S,n — 2, 


On = 


On = 


(5.21) 


; j = l,2,...,n-4. 


( 1 ^ 


n{n + 2 ) 

-i( 2 n 


n{n + 2 ) 


; j = n-l,n. 


Note that is defined in (2.10) determine the Frobenius expansion around the 
boundary for the equation (5.18). The thermalizing Green function can 

be read out by by applying the matching condition, (2.7), on This is because Hy 
can be obtained by integrating py. This just introduces a power of z which adjusts 
to give the right A+ root for the fluctuation Hy. Now Hy is related to the vector 
fluctuation hy^i which reduces to at the boundary v equals the boundary time 
t. Therefore reading out the time dependence of py is sufficient to extract the time 
dependence of the the two point function {TtxiTt^i). 


Green function and quasinormal modes 

Using the equations (5.20),(5.21),(2.14) and (D.4), with u = 1, the moments Mi 
are obtained. Then the retarded two point function of (u,a;^) component of the 
stress energy tensor of the boundary theory, is evaluated using (2.18). We used 
Mathematica to evaluate up to 3500 moments ^ 2 ( 0 ;), the results converged at this 
order of moments. The Green function is constructed to (9(u^®°°) in future time. 
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Figure 7: The red dots show the value of the hydrodynamic frequency obtained 
from the slope of the long time behaviour of the Green function (at ti = —5) for 
various values of k. The blue curve plots the expected behaviour, —Im u = \k'^. 

V and then we approximated it by its (1750|1750) Fade approximant. From these 
moments the Green function is evaluated for various values of the momentum k. 

Figure 6a shows the thermalizing Green function (solid blue curve), which starts 
close to the vacuum Green function (dot-dashed red curve), and deviates away from 
it within one horizon time. 

The hydrodynamic pole of the thermal two point function or the quasi¬ 

normal mode of the vector perturbations of the metric obeys the dispersion relation 

uj = -i^k^. (5.22) 

rj, s and T are the shear viscosity, entropy density and temperature, respectively. 
The universal value of rj/s for the vector mode of metric perturbation is known to 
be (dvr)”^. The temperature of the AdS^ black hole is T = 7r“^, in which the radius 
of AdS and that of the horizon are normalized to unity. Therefore the dispersion 
relation for the hydrodynamic quasinormal frequencies is given by 

w = -i — . (5.23) 

As discussed in section 4 the long time behaviour of the thermalizing Green function 
is dictated by the lowest quasinormal mode. Therefore it is possible to extract the 
value of the —Imcu from the slope of the logarithmic plot of the thermalizing Green 
function for k « 1 and show that it obeys the dispersion relation given in (5.23). 

In figure 6b, the logarithmic plot of the absolute value of the thermalizing Green 
function (solid blue curve) is plotted along with the line measuring the slope of its 
large time behaviour (dot-dashed red line). The slope of this straight line gives the 
value of —Im uj. The plot is for fixed values of ti = —15 and k = 10. Figure 7, shows 
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the red dots for values of —Im u for several values of k obtained from Mathematica. 
For small fc, these values lie very close to the expected parabolic behaviour (solid 
blue curve), predicted by the hydrodynamic dispersion relation (5.23). For larger 
k, these values lie inside the parabola. This behaviour agrees with figure 14 of [27] 
which studied this quasinormal mode for large values of k in the AdS^ black hole. 

6 Conclusions 

We have developed a recursive method to obtain the time dependent Green functions 
in the thin shell AdS Vaidya background. Using the intuition developed from this 
method we showed that the long time behaviour of the Green function is determined 
by the lowest quasinormal mode of the corresponding black hole. Thus our analysis 
provides an explicit realization of the general conjecture made in [34] that time scales 
in thermalization are determined by the quasinormal modes. We applied the method 
to study Green functions in thin shell Vaidya geometries in AdS^ and AdS^. Using 
this method we obtained the universal ratio of shear viscosity to entropy density by 
studying the relaxation of the time dependent Green function of the vector metric 
perturbation in the AdS^ Vaidya shell. 

The methods developed in this paper to study time dependent Green functions 
can be generalized for fermion Green function which was obtained by direct numerical 
integration in [24]. It will be interesting to generalize this method to Vaidya shell 
geometries which have a finite width found in [13]. These geometries are a more 
accurate description of the thermalization process. The dependence of the width 
of the shell on the time dependence of the Green function will be interesting to 
extract. It is interesting to study thermalization of Green function of higher spin 
fields to study the spin dependence in thermalization. It is known that a non local 
probe like entanglement entropy is the slowest to relax to equilibrium [21]. It is will 
be interesting to see if higher spin fields relax faster or slower compared to scalars 
and entanglement entropy. In this context the exact solutions to wave functions of 
higher spin fields found in the BTZ background [38, 39] will prove to be useful to 
obtain analytical results. Finally it will be useful to use these lessons learned in 
holography to understand time dependent Green function in conformal field theories 
during thermalization along the lines developed in [40] . 
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A Scalar wave functions in BTZ 


In this appendix we derive the scalar wave functions in the BTZ geometry as an 
expansion around the boundary by transforming the known wave functions which 
are written as expansions around the horizon. Consider the wave functions given in 
(3.12) which admit a natural expansion near the horizon. We can use the following 
transformation property of the Hypergeometric function to write it as an expansion 
around the boundary, 


F(a, 6, c, z) 


P(c)P(c — a — b) 
Y(c-a)P(c-6) 


F(a, h^a + h — c + 1,1 — z) 


(A.l) 




,-a-b '^{c)T{a + b - c) 
T{a)T{b) 


F{c — a, c — b, c — a — b + 1,1 — z) 


Substituting (A.l) in (3.12) we obtain 
(;2, V, 0) = X 

Bi{iy,u,k)z^-F Q(A 

+ B2{u,u,k)z^+F 


- - i{u - k)), + k)),/^-,z 

- i{u + k)), ^(A+ - i{u - k)), A+, 


(A.2) 



=e 


_ —iLOV+ik^) 


(1 - zf^^ X 

A_ 


(A.3) 


B^{iz,uj, k)z^-F ( ^(A_ + i{uj + k)), ^(A_ + i{uj - k)),^-,z‘^ 


+ Bi{v,u,k)z^+F ( ^(A+ + i(a; - k)), ^(A+ + i(a; + k)),^+,z^ 


The coefficients Rj, i = 1,... 4, are made up of gamma functions obtained from 
(A.l), and do not depend on 2 ;. The hypergeometric functions of (A.3) can further 
be related to the hypergeometric functions of (A.2), by the transformation, 

F(a, b, c, z) = {1 — zy~°‘~^F{c — a,c — b,c, z). (A.4) 
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Thus, the only two independent near boundary solutions are, 


\z,v,4>) 




F Q(^- - - ^))>+ k)), A_, 2;^ ) , 

_ ^-iLLjv+ik(l> 


(A.5) 


^ ( ^(A+ - i{uj + k)), ^(A+ - i{u - /c)), A+, 2 ;^ 


B Thermal Green function in BTZ 

The thermal Green function is obtained by taking the Fourier transform with respect 
to oj of the following thermal Green function in Fourier space [23] , 


^thermal 

R 


{k,u) 


r ( ^(A+ + i{uj + k))\ F Q(A+ + iiu - k)) 


2 sin(7rA+)(F(A+))2 
X (cos(7rA+) cosh(7ra;) — cosh(7rA;) — isin(7rA+) sinh(7ra;)) . 


(B.l) 


The Fourier transform is taken by performing integration over and 

completing the contour in the lower half cn-plane. It can be shown that the integral 
vanishes over the arc of sufficiently large radius. The result is a sum over the residues, 
which are evaluated at the poles of the gamma functions in the lower half ca-plane. 




2 sin(7rA+)(F(A+))2 


(B.2) 


X 


E E kAe--;(«-<.)rQ(A++ *(!.,»-„«:))[) 

n a=±l ' 


X r Q(A+ + + ak))^ V Q(A+ - i(a;“ + ak)) 

X [cos(7rA+) cosh(7ra;“) — cosh(7rfc) — i sin(7rA+) sinh(7ra;“)]. 
These poles, a;(n, a), are the quasinormal modes of the BTZ black hole, 

a;“ = afc — i(A++ 2n), a = ±1. (B.3) 


C Mixed Fourier transform of Green function in AdS 

In this appendix we solve equations of the form (5.18) with the analog of the boundary 
conditions in (2.5) hrst in the frequency-momentum space. We then perform the 
partial Fourier transform in the frequency space. This results in the Green function 
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in AdS which satishes the boundary condition of (2.5). Consider (5.18), with a more 
general mass term, 


z‘^p” + z{—3 + 2iujz)p'^ — (—3 + + k'^z'^ + 3iojz)p^ = 0, (C.l) 

this equation has the following Fourier space solution, where v = a/1 + 

p^{k,u,z) = e~^‘^^z'^ + B{u})J_iy (^zVoj'^ — • (C^-2) 

Using the property of the Bessel function for small argument, the following near 
boundary behaviour of the solution is obtained. 


// A ^^0. 2/ 4^ 4 fzVuJ^ -k^\ 1 , f Zy/uF^W\ " 1 

Py{k,u,z) -)■ 2 ; 74 ( 0 ;) --- —- r + B{u) 


r(i + 0 


r(i - 0 

(C.3) 


The delta function boundary condition in (2.5) in the time domain requires that near 
the boundary the wave function in Fourier space solution goes to 


Py{k,u,z) 


z-s-O J- A- 

- > — • ^ 
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This determines B{cj) to be 


B{u7) 


F(l-n) 

^ V 2 ) 


(C.4) 


(C.5) 


where, A± = 2 ± u. 

We now use the following asymptotic behaviour of Bessel functions to obtain the 
behaviour of the solution near the origin of AdS, 


jV ( 2 ^) 
>4—1/ ( 3 ^) 



2/1^71 

COS [x - 

TTX \ 2 

cos(7rz/) cos 

TlX L 





— sin(7rz/) sin 



(C.6) 



As z —)■ 00 , the solution. 


p^{k,u,z) 




(A(a;) + B{u) cos(7ri/)) cos -^ 

/ TIU 7l\' 

- B{u) sm(7ri/) sm - — - - j J, 


71 

4 


(C.7) 


where, x = z\/u^ — k^. The solution near the origin has both ingoing and outgoing 
behaviour due to the presence of cosine and sine. In order to impose ingoing boundary 
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condition at origin, the solution is rewritten as a linear combination of ingoing and 
outgoing solutions, 


where C and D in terms of A and B, are 

C(u,) = t (/IH + B(u,)e -‘"), (C.9) 

D(aj) ^ AaIoj) + B(u)e‘"'). (C.IO) 

The ingoing boundary condition at origin, means that B>(u) = 0, hence, 

A(u) = (C.ll) 

Substituting this in (C.2), we obtain, 

cn, z) = B{u)e-^^^z^ [-e^^’^J^{x) + J-u{x)] (C.12) 

= -ism{nu)e^^''B{uj)e-^‘^^zH^^l{x). 

The Fourier transform of the above equation with respect to u, gives 

GPHv,k,z)= (C.13) 




using the integral representation of H S(^), 


H^-l{x) = -e 

TT 


2 X 


dy exp 


1 X 

■X h/ + — 

2 V y 


y^-\ 


(C.14) 


Gf\v,k,z) = 


ITTl/ 

e 2 


—z 


2'^27rr(z/ 

/ OO 

du exp 

■OO 

_• TT inu 

e 4e 2 

2'^27iT{h'' 


^2-u 


u^z^ 


22 / 


dyy'" ^exp 


+ u{v + z) 


z'^k"^ 


(C.15) 




_ • ITTl/ 

e 40 2 


2^27rr(i/ 

1 

2^r{u)^r 


V^z ^-’'2 


dy y'^ 2 exp 
kz^ " 


■ /^2n fv 


y/v^ + 2vz 


+ 




y 


K 


— I' 


-iky/v"^ + 2vz) 


y 2 +iy 


k 


+ 2vz 


u+i, 


H^_^l_i{ky/v'^ + 2vz). 
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In the last step, K\{x) = has been used. The retarded two point 

function is obtained from the Feynman two point function using the relation 

GRit, t') = 9{t - t') {Gpit, t') + t')). (C.16) 

The following properties of the Bessel functions, 

h‘^\xY =Hf{x), (C.17) 

= icsc(7rA)(e-*"^ - e*"^)JA(a:), 

are used in last two equations, to obtain the retarded Green function which is given 
by 

Gf%v -h,k,z)= C0(v - ( ■ *; ==) (C.18) 

\^{v - tiY + 2{v - ti)z J 

X {v - tiY + 2{v - ti)z^ , 

where, G = Though we started out by assuming u is not an integer and the 

two independent solutions are given by (C.2), the discussion can be generalized to 
the case when u is an integer, leading to the same hnal result given in (C.18). A 
simple way to see this is that in the hnal result the order of the Bessel function is 
fractional. 

D Recursion relation for retarded Green functions in AdS 

In this appendix for completeness we obtain the recursion relation for the mixed 
Fourier space retarded Green function in AdSd+i at n = 0. This is useful to imple¬ 
ment the recursive algorithim to obtain the thermalizing Green function in thin shell 
Vaidya geometries. The solution for a scalar held in mixed Fourier space in AdS^+i, 
has the following dependence on the Bessel function, 

= 0,k,z-,ti) = z^+R{z). (D.l) 

Here we have assumed H < 0. 

where, from (3.3) and (C.18), we see that, A+ = d/2 + v. The diherential equation 
satished by R{z) is, 

(ti - 2z)R!\z) + 2(A - l)R!{z) + kHiR{z) = 0. (D.3) 
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Substituting the ansatz, R{z) = in the above equation, the following 

recursion relation for is obtained, 


_ 2(?7- 1)(?7- 1 X)Jn-i k‘^tiJn-2 _ for 71 > 9 

tiuyn — 1 ) 


(D.4) 


To get the solution of (D.2) from the above recursion relation, Jq and Ji are chosen 
to be, 

k 

\ti\ 
k 
|ti 

Thus we obtain the Frobenius series expansion of AdS scalar field solution, around 

z = 0, 


^+1 



(D.5) 



kJ_^_3{\k\\ti\). 

(D. 6 ) 


= 0, /c, 2 ;; ti) = ^ JnZ^' 

n=0 


(D.7) 


E Details of Mathematica files 

The recursion method of section 2 is used to numerically construct the Green function 
using Mathematica for the three cases considered in this paper. The hies for each 
case is attached. 

1. ads3-green.nb: The Green function for scalar held in AdSs Vaidya is numeri¬ 
cally constructed to order 56 in moments and hgures 2 and 3 are obtained. 

2. shear-green.nb: The Green function for shear metric perturbations h^ix^ 7 ^ 0 
in AdS^ Vaidya is numerically constructed to order 200 in moments and hgure 

5 is obtained. 

3. vector-green.nb: The Green function for metric perturbations hxix^ 7 ^ 0 in 
AdSs Vaidya is numerically constructed to order 3500 in moments and hgure 

6 is obtained. 
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